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Abstract : In the absence of accretion torques, a pulsar in a binary syctem will 
spin down due to electromagnetic dipole radiation and the spin-down power wi^ 
drive a wind of relativistic electron-positron pairs. Hinds from pulsars wit£short 
periods tai 1 1 prevent any subsequent accretion but may oe confined by the corndKion 
star atmosphere, wind or magnetosphere to fora a standing shock. The authors Hk 
investigate the possibility of particle acceleration at such a pulsar wind 
and the production of VHE and UHE gamma rays from interactions of accelerate^^Bt 
protons in the companion star's wind or atmosphere. They find that in close ^Hg 
binaries containing active pulsars, protons will be shock accelerated to a mwa 
energy dependent on the pulsar spin-down luminosity. If a significant fraction u. 
the spin-down power goes into particle acceleration, these systems should be 
sources of VHE and possibly UHE gamma rays. The authors discuss the application of 
the pulsar wind model to binary sources such as Cygnus X-3, as well as the 
possibility of observing VHE gasna-rays from known binary radio pulsar systems. 
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Abstract 


In the absence of accretion torques, a pulsar in a binary system will spin down due to 
electromagnetic dipole radiation and the spin-down power will drive a wind of relativistic 
electron-positron pairs. Winds from pulsars with short periods will prevent any subsequent 
accretion but may be confined by the companion star atmosphere, wind or magnetosphere 
to form a standing shock. We investigate the possibility of particle acceleration at such 
a pulsar wind shock and the production of VHE and UHE gamma rays from interactions 
of accelerated protons in the *ompanion star’s wind or atmosphere. We find that in close 
binaries containing active pulsars, protons will be shock accelerated to a maximum energy 
dependent on the pulsar spin-down luminosity. If a significant fraction of the spin-down 
power goes into particle acceleration, these systems should be sources of VHE and possibly 
UHE gamma rays. We discuss the application of the pulsar wind model to binary sources 
such as Cygnus X-3, as well as the possibility of observing VHE 7 -rays from known binary 
radio pulsar systems. 


I. Introduction 


Magnetized, non-accreting neutron stars will spin down through energy losses due primarily 
to magnetic dipole radiation (Ostriker and Gunn 1969). Since the observed radiation from 
pulsars at radio frequencies and above accounts for only a small fraction of their total spin 
down energy loss, it was originally proposed that the bulk of this energy is carried away in the 
form of large amplitude electromagnetic waves at the rotation frequency of the pulsar. The 
development of pulsar models since the work of Ostriker and Gunn has shown that vacuum 
conditions outside the star cannot exist in the presence of electric forces many times the 
gravitational force at the stellar surface (Goldreich and Julian 1969). Furthermore, it is very 
likely that pair production takes place somewhere in the pulsar magnetosphere (Sturrock 
1971), raising the particle density in these regions to several orders of magnitude above the 
corotation charge density. Under these conditions, the plasma frequency at the light cylinder 
will be significantly higher than the frequency of the* dipole radiation (see e.g. Arons 1981). 
Consequently, the dipole electromagnetic waves cannot propagate and the pulsar spin-down 
energy will most likely be transported by a relativistic MHD wind consisting mostly of 
electron-positron pairs (Rees and Gunn 1974). 

Confinement of pulsar winds can occur if there exists enough pressure surrounding the 
pulsar to balance the wind ram pressure. The static pressure of the interstellar medium 
is not large enough to confine the wind of an isolated pulsar. However, if the pulsar is 
surrounded by a supernova remnant or is moving relative to the interstellar medium, then 
the wind can be confined. In their model for the Crab nebula, Rees and Gunn (1974) 
proposed that the expanding shell would confine the pulsar wind at a standing shock where 
the wind ram pressure balances the total magnetic field and particle pressure in the nebula. 
They assumed that since the field and particles in the nebula have come from the pulsar 
(the present field in the Crab nebula is many times larger than the remnant field of the 
expanded stellar envelope), the nebular pressure results from the energy accumulated over 
the pulsar’s lifetime. Kennel and Coroniti ( 1984a, b) have extended the ideas of Rees and 
Gunn to construct a detailed pulsar wind model of the Crab nebula. Cheng (1983) has 
considered the confinement of pulsar winds by ram pressure from the pulsar’s motion through 
the interstellar medium. In this case a bow shock forms ahead of the star at the pressure 
balance point and the synchrotron radiation from the shocked wind particles is predicted to 
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be observable in X-rays. A bow shock of this kind may have been observed in Ha emission 
surrounding PSR1957+20 (Kulkarni and Hester 1989)- 

In this paper, we consider the possibility of confinement of a pulsar wind by a binary 
companion. Pulsar winds in binaries were investigated by Davies and Pringle (1979) and 
were first proposed as a model for Cygnus X-3 by Bignami et al. (1973). More recently, the 
interaction of a pulsar wind with a binary companion has been conside. “d in models for the 
evolution of low-mass X-ray binaries (Ruderman, Shaham and Tavani 1989, Ruderman et 
al. 1989), in particular for the eclipsing 1.6 ms pulsar PSR1957+20 (Phinney et al. 1988, 
Kluzniak et al. 1988, Cheng 1989). In Section II, we discuss the conditions under which 
pulsar winds can form in binary systems and the confinement of the pulsar wind either by 
gas or magnetic pressure in the companion star atmosphere or by the ram pressure of the 
companion star wind. In both cases, the location of the stationary shock front relative to 
the separation of the two stars is calculated. In Section III, we explore the possibility of 
first order Fermi acceleration of charged particles by the pulsar wind shock and estimate the 
energy to which the particles can be accelerated. The high energy gamma rays resulting from 
interactions of these particles with material in the companion star atmosphere is discussed 
in Section IV. In Section V we discuss the application of these results to binary sources such 
as Cygnus X-3 which have been observed in UHE and VHE gamma rays, as well as the 
possibility of observing this type of emission from known binary pulsar systems. 


II. Formation of the Pulsar Wind Shock 


The power in magnetic dipole radiation from a pulsar with rotation frequency f! and magnetic 
dipole moment m is 




2m 2 fi 4 sin 2 9 
3c 


4 x 10 43 erg s _1 B 2 12 P^ 


( 1 ) 


where is the period in ms, B\-„ = (f? o /10 12 Gauss) is the surface magnetic field and 9 is 
the angle between the dipole and rotation axes. We assume that all of this power appears as 
a relativistic wind which carries both particles (predominantly electron-positron pairs) and 
wound-up magnetic field away from the pulsar. Since the magnetic field is dipolar (as r~ 3 ) 
inside the pulsar light cylinder, r^c = c/tt = 5 x I0 6 cmP mi , and toroidal (as r _1 ) in the 
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wind, the field strength at a distance r is 

B = = 8 x lO’Gauss BuP^f-f), (2) 

ric r r 

where r 0 = 10 6 cm is the neutron star radius. Eqn (2) implicitly assumes equipartition 
between magnetic field energy and particle energy in the wind (Auriemma, Gaisser and 
Lipari 1988). Winds from pulsars with short periods have higher magnetic fields because 
the light cylinder, inside which the field falls off most rapidly, is closer to the neutron star. 

Accretion onto the neutron star will not take place as long as the pulsar light cylinder, tlc , 
where the wind would be generated, is inside the Alfven radius, = 1.5 x 10 8 cm By? Mi^ 1 . 
This will be the case for pulsars with periods 

P < 31 ms Bll 7 M~ s 2/7 , (3) 

where M 18 = Mj 10 18 gs _1 is the accretion rate (see also Ruderman, Shaham and Tavani 
1989). In these cases where the light cylinder is inside the Alfven radius, the ram pressure 
from the wind, Lj/47rr 2 c, everywhere exceeds the ram pressure of the accretion flow, pv 2 . 
In order to have a stable force balance between the pulsar wind and the accretion flow, pv 2 
must fall off vith r more slowly than the 1 /r 2 dependence of the pulsar wind. For spherical 
accretion, pv 2 sa 1/r 5 ^ 2 , while tor disk accretion, both in a disk (Shakura and Sunyaev 1973) 
and thick disk models pv 2 falls off faster than 1/r 2 . Hence, there can be no confinement of 
the pulsar wind by an accretion flow. In fact, if the pulsar light cylinder is inside the Alfven 
radius, there can be no accretion flow at all since the pulsar wind is capable of blowing away 
accreting material at all radii (except where the companion star dominates the gravitational 
potential). 

In the absence of an accretion flow, there are several other possibilities for confining 
the pulsar wind in a binary system, at least over a limited solid angle. We consider three 
general cases: 1) confinement by the static atmosphere of the companion where the gas 
pressure balances the ram pressure of the pulsar wind 2) confinement by the companion star 
wind where ram pressure of the two winds balance 3) confinement by the companion star 
magnetosphere. 

In the second case, the wind from the companion star may be driven either by processes 
intrinsic to the companion or by the relativistic wind from neutron star (or by a combination 
of the two). Figure 1 shows a schematic representation of the pulsar wind shock formation. 
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A reverse shock will form on the pulsar side of the contact discontinuity which separates the 
pulsar wind from the companion star wind (or atmosphere) at the pressure balance point. We 
propose that particles may be accelerated at this shock. A reverse shock in the companion 
wind (as shown in Fig. 1 of Phinney et al. 1988) may also form if the wind speed becomes 
supersonic, but as pointed out by Chen" (1989), need not exist in all cases. The extent to 
which field line reconnection occurs in the shocked pulsar wind will probably depend on the 
geometry of the pulsar wind flow around the companion and could be an additional source 
of accelerated particles. 


a) Companion Star Atmosphere 


In the case of a static atmosphere, Ld/(4nr 2 c) is balanced by the gas pressure, pkT/mn, 
close to the companion star so that the shock distance from the pulsar, r = r, ss a — r.. In 
an isothermal atmosphere, the density has an exponential dependence, p = p 0 exp(-z/h ), 
so that significant changes in density occur over distances comparable to the scale height 
h = kT I mug. To determine the density near the shock, we need to know the temperature in 
the atmosphere. A rough estimate may be obtained from requiring that the temperature be 
high enough to radiate away the pulsar wind power absorbed by the star, L, = eLj, where 




(4) 


and fl. is the solid angle subtended by the companion. This condition, S.crT 4 = eL * , where 
5, = 27rr 2 (l — r„/a) is the area oi the companion which can absorb wind power, gives the 
following temperature and scale height: 

£ l/4 
r. 


T = 2.2 x 10 6 KB 1 1 2 /2 / j -, 1 (^®) 1/2 


h = 5.9 x 10 9 cm B \? 


(1 — r./o) 1 / 4 
c l/4 


( 5 ) 


( 6 ) 


R q ' y M '(l-r./a)‘/4 
where a is the sum of semi-major axes of the binary system (i.e. the separation of the stars), 
and M and r. are the companion star mass and radius. The atmospheric density at the 
discontinuity is then determined by the force balance condition to be 


,. = 1.2 x lO-Vm- 3 ^,^) 3 '’ ' 


- 1/4 


(r./a) 2 


r„ (l-r./a) 7 / 4 ’ 


( 7 ) 
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b) Companion Star Wind 


In the case where the pulsar wind is confined by a stellar wind from the companion, we 
balance Li/(Airr 2 c) with (pu 2 ) w = M m Vy,/(4irr 2 ). The shock radius is thus dcteimined by 


L d 

Citric 


M-atV 0 


1 - 


r. 


a - r,J 


1/2 


( 8 ) 


4?r(a — r,) 2 

where A/ w and v x are the mass loss nte and terminal velocity of the wind, and r, is the 
distance of the shock from the neutron star (cf. Fig. I). Note that Eqn (8) is valid only 
for shock radii well outside the critical point radius, where the wind speed equals the local 
sound speed. 

From Equation (8), (pu 2 ) w reaches a maximum value at r^x = a— |r,. The condition that 
{pv 2 ) w balance Lhe pulsar wind somewhere between the two stars, and thus the requirement 
for pressure balance, is {pv 2 ) mtiX > Li/Axr^c or from Eq. (8) 

5\ 2 , r„ • , Uoo x , M 




( 9 ) 


where u« c = (2 GM/r,) 1 ? 2 is the escape velocity from the companion. If Eqn (9) is satisfied 
then Eqn (8j gives a solution to the shock radius r, and we require that r^c < r s < a — r,. 
Figure 2 shows the location of the shock from numerical solution of equation (8) as a function 
of r,/a. In the limit r g << a — r,, the approximate analytic solution is 

r , _ \/K 


4^ — - 


L d 


a 1 + y/x: 
x10 7 £j 2 


2x10 ’’B&Pzi l-fc)'/* 


(10) 


(U) 


Af w v w c 

The quantity A w is essentially the ratio of pulsar wind ram pressure to stellar wind ram 
pressure and is therefore the critical parameter determining solutions for r,. 

From mass flux conservation, the stellar wind density near the discontinuity is 

P. = 2-64 x 10' ,3 cm' 3 r 22 52 , (12) 


i r*/a 

‘- 0 = 5*75 


(l-r,/a) J 


where we have used the approximate solution of Eq. (10). 

In the case where the stellar wind is driven by the companion star’s self-generated ra- 
diation pressure, M w and are free parameters. In the case where the 6tellar wind is 
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induced by absorption of pulsar wind luminosity, solutions are almost completely defined by 
the pulsar wind power, which determines the ram pressure of both winds. We consider two 
special cases of pulsar-induced winds. 


!. Radiation Pressure Driven Induced Winds 

The pulsar wind power absorbed by the companion star, Z. « eLj is greater than the 
Eddington limit, Z-Edd = 4nGcM/K, when 

(^) (0 > 3.15 x 10-«, (13) 

where k it the opacity (with k c = 0.4 cm 2 /g the electron scattering opacity) and e is defined 
in Eqn (4). 

If tLd > Leu . then heating by the pulsar wind can drive a wind from the companion 
through radiation pressure. The wind terminal velocity will be 

v oo — [2GA/(r — 1 )/r,] 1 ^ 2 , (14) 


where T = tLiKjMcGMc is the ratio of the radiation pressure force to the gravitational force. 
The induced mass loss can be estimated by the condition = eLj (assuming that all 

of the absorbed energy goes into bulk kinetic energy of the wind with terminal velocity v,*,). 
In the limit T » 1, 


Mi = 


47TCT. 

K 


= 6.6 x 10 ?2 gs -1 



(15) 


and M\ is independent of Lj. If the opacity driving the wind comes mostly from lines, as 
is often the case in stellar atmospheres (Mihalis 1978, Castor et al. 1975), then k may be 
much higher than the electron scattering opacity. Effects of seif-shielding by the stellar wind 
may also reduce M\ (Eichler and Ko 1988). Equations (14) and (15) give 


Mis (^) = 3.7 x 10 7 P- 2 5, 



(16) 


Figure 3 shows numerical solutions of Eqn (8) for the shock location r,, using Eqn (16) in 
the expression (II) for A w . In the limit r t « a — r., the variation of r,/a in Fig. 3 comes 
from the e dependence, which reflects the relative amount of pulsar wind power absorbed by 
the companion. The solutions shown in Fig. 3 must also satisfy eLj > Z-Edd- Some M = Mi 
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solutions may be unstable in the sense that Mi could be large enough to turn off the pulsar 
wind if the Alfven radius of the induced mass accretion onto the neutron star is inside the 
pulsar light cylinder, violating Eq (3). It is difficult to give a simple condition for systems 
with stable solutions, since the relation between the accretion rate M and the stellar mass 
loss rate A/ w depends on the geometrv of the system, wind angular momentum, etc. 


2. Thermally Driven Winds 

If eLd < L e<w, the pulsar may still induce a wind from the companion star if the pulsar 
wind can efficiently heat the atmosphere to a temperature at which the mean particle ve- 
locity exceeds escape velocity. Kinetic energy gained by particles in the atmosphere can be 
converted directly to wind kinetic energy if the expansion time is less than the radiation 
(bremsstrahlung) timescale. Thermal wind models have been studied in th~ context of ac- 
creting sources where thermal X-rays from the neutron star drive a self-sustaining wind from 
the companion star. London, McCray and Auer (1981) and ^London and Flannery (1982) 
have shown that the resulting mass loss from X-ray heating is approximately 

M ~ 5 x 10 _,, g s _I <^“ 1/2 L X , (17) 


where (j> = GM/r, = v eK f 2 is the surface gravitational potential and L x is the X-ray lumi- 
nosity incident on the companion. Substituting this expression into Eqn (11) and assuming 
v w — v 'sc, the parameter A w which determines the shock radius becomes 

0.466 


At a — 


« ©' 


(18) 


where we have taken L x = (Lj- 

Thermally driven winds may also occur in some of the radio pulsar binary systems. Cheng 
(1989) has investigated a self-consistent model for pulsar-induced stellar winds, in which the 
estimated mass loss rate from the companion is given (from his Eqn [17]j by, 

- 1/2 




( - 


(19) 


where P. is the pulsar wind ram pressure, A 0 is the cross sectional area at the base of the 
wind flow (?» irrl), ( = (u oc ,/u egc ) 2 and fi, 7 c and 6 are dimensionless parameters specifying 
the flow geometry. From the above expression (with fi = ] .5, = 5/3 ,6 = 0.5,), we can 
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derive the product of mass-loss rate and terminal velocity, 

£) - 1.7 x £)* (£)’ /J £)’. (20) 

Inserting the above expression into Eqn (11) gives a particularly simple form for the param- 
eter A w which determines the shock location: 

3 / „ \ 2 


A w = 11.5 


(=) ©' 


( 21 ) 


Solutions for r t as a function of r ,fa are therefore dependent only on the wind terminal 
velocity. Figure 4 shows numerical solutions of Eqn (8) for the shock location in the ther- 
mal wind model of Cheng (1989) as a function of r./a and (voe/ u e*c)- The solutions are 
constrained by Eqn (9) to require 

Sr.\ ,p 


in the limit that the shock radius lies well outside the critical point. 


( 22 ) 


c) Companion Star Magnetosphere 


If the companion star has a surface magnetic field, B m , then magnetic pressure may be 
sufficient to stand off the pulsar wind. In this case, the shock location can be found by 
balancing [B„(r*)] 2 /8 , r with L d /(Axr*c). Assuming a dipole field, r, is determined by 


Bl f r. y = L d 
8tt \a — r s / 4 7rr 2 c 


(23) 


Solutions for the shock location from Eqn (23) are determined by the quantity 

2 L d 

m BWc 


(24) 


which scales with the ratio of pulsar wind ram pressure to magnetic pressure. Figure 5 
shows numerical solutions for r, from Eqn (23) as a function of the parameter A m . The 
condition that the shock form above the companion star surface, or ( a — r s ) > r., requires 
that A m < (1 — ?“./a) 2 < 1, which from Eqn (24) gives a limit on the companion star field 
capable of standing off the pulsar wind, 


B.> 



i___ 

(a - r.V 


(25) 
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III. Shock Acceleration of Particles 


Once a shock forms in the nulsar wind as a result of confinement, the conditions may 
be favorable f A r first-order Fermi acceleration of particles. Shock acceleration has gained 
considerable attention recently as a mechanism for generating highly energetic particles in 
a variety of sources which include supernovae (Lagage and Cesarsky 1983), the interstellar 
medium (Blandford and Ostriker 1978), the Earth's bow shock (Ellison and Eichler 1984), 
stellar winds from hot stars (Cesarsky and Montmerle 1983) and pulsar winds in supernova 
remnants (Gaisser, Harding and Stanev 1987, 1989). If the shock is formed by collisionless 
processes, particles can travel back and forth across the shock front by scattering from 
magnetic irregularities, gaining some energy cn each crossing. The spectrum of accelerated 
particles escaping downstream is a power law with index dependent on the shock compression 
ratio. The theory of diffusive shock acceleration (cf. Drury 1983 for review) has concentrated 
primarily on strong parallel shocks (where the magnetic field is parallel to the shock normal), 
in which the jump conditions lead to a compression ratio 



4 

1 + 3/M 2 


4. 


(26) 


Here M is the Mach number of the shock. The standard treatment also is for non-relativistic 
shocks for which tii, u 2 c. Such a shock produces a power law spectrum with index of 
the differential energy spectrum a = 2 + 4 jM 7 . 

The shock in the relativistic pulsar wind differs from tht canonical case in two ways. 
First, because of the toroidal field, the shock is quasi-perpendicular rather than parallel, so 
that the jump conditions require compression of the magnetic field across the shock. This 
actually tends to increase the acceleration rate because the particles spend more time in the 
vicinity of the shock and gain energy by drifting through the potential drop induced along 
the shock (Jokipii, 1987). Second, the pulsar wind shock is relativistic in that the velocity 
of the unshocked wind relative to the shock front is u t ts c. Treatments of acceleration 
by relativistic shocks (Peacock 1981; Kirk and Schneider 1987) indicate that they are more 
efficient at accelerating particles than non-relativistic shocks (the energy gain per crossing is 
larger). Recent Monte-Carlo simulations of Ellison, Jones and Reynolds (1990) give a shorter 
acceleration time in relativistic shocks and show that the accelerated particle spectrum is 
flatter than for a non-relativistic shock with the same compression ratio. However, the 
resulting spectral index is somewhat uncertain, since the two approximations usually made to 
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treat scatterings in the non-relativistic case, pitch angle scattering or hard-sphere scattering, 
give different results in the relativistic case. 

For the present, we simply assume that this configuration of a relativistic shock with 
a quasi-perpendicular field configuration is capable of accelerating particles at least as ef- 
ficiently as the non-relativistic, parallel shock. We then use the well developed theory of 
non-relativistic, parallel shocks to estimate the acceleration rate and hence the maximum 
energy accessible. The acceleration rate for a relativistic particle is (Lagage and Cesarsky 
1983) 


3f(f+l) 0 ’ 


(27) 


where D is the diffusion coefficient (which we assume to be the same both upstream and 
downstream) and U\ is the flow velocity of the upstream (unshocked) fluid into the shock. 
We now use the minimum value of the diffusion coefficient to obtain an upper limit on the 
acceleration rate and hence an upper limit on the maximum energy that can be achieved. 
We take Anin = r^v/Z, where n, = 3.3 x 10 9 cm£x e v/l? is the Larmor radius, v is the 
particle velocity, and E is the particle energy. For Uj fts c % v the resulting estimate of the 
acceleration time is 


t a a 


£ €K + i) 


(28) 


ceB ({-1) ‘ 

Because of the inverse dependence on field strength, which from Eq (2) is large, and the 
relativistic velocity of the wind, the pulsar wind shock is extremely efficient at accelerating 
particles to high energy. 


a) Conditions for Shock Acceleration 

The conditions under which diffusive shock acceleration can take place are 1) the plasma 
must have 0 < 1, where fi = Blfin^pc 7 is the ratio of magnetic to particle energy density 
in the plasma. This condition is necessary for the existence of a strong shock (Kennel and 
Coroniti 1984a) and also for the development of magnetic turbulence which scatters the 
particles 2) the shock must be collisionless, forming as a result of magnetic turbulence or 
other collective plasma effects rather than by collisions of individual particles, i.e. the plasma 
and photon densities both upstream and downstream of the shock must be low enough to 
prevent inelastic collisions from dominating the scattering. 

The first condition is equivalent to requiring that the Alfven velocity be less than the 
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particle velocity, which for relativistic flows is close to c. Kennel and Coroniti (1984a) and 
Enunering and Chevalier (1987) have examined relativistic pulsar wind solutions as a model 
for the Crab nebula. In these models the solutions are constrained by matching the jump 
conditions at the shock 'and by the expansion velocity of the shell which confines the wind. 
The free parameter, o, characterizing the flow solutions in these models is the ratio of the 
magnetic energy flux to particle energy flux in the unshocked wind and thus is the same 
as the parameter defined above. The value of a (or 0) in a pulsar wind is not well 
determined, because there exists no global solution which relates the particle production 
in the pulsar magnetosphere and flow of these charged particles through the light cylinder 
with the wind solutions outside. Large values of a, where the flow is highly magnetized, 
produce weak shocks with low thermalization, while small values of a give a strong MHD 
shock approaching the classical hydrodynamic limit. The above studies find that a very 
small value of a (ss 10" 3 ) is required to match the present expansion rate and synchrotron 
luminosity of the Crab nebula. There are indications therefore that pulsar winds have $ < 1, 
favorable for shock acceleration. If particle energy flux dominates magnetic energy flux in 
a pulsar wind, then the strength of the magnetic field at the shock would be lower than 
that derived by assuming equipartition (Eqn [2]) by a factor of at, reducing the acceleration 
energy by the same factor. 

The second condition requires that the mean free path for inelastic scattering or absorp- 
tion processes be large compared to the acceleration length, l a = ct a . To examine whether 
the second condition is satisfied in the pulsar wind, we will consider both nuclear collisions 
and photo-pion production in the radiation field of the companion. The nuclear collision 
mean free path is /„, sb Ap P fp, where App a 60 gem -2 is the interaction length in H at 
around 1 TeV (App decreases slowly with energy to a value of 35 gem" 2 at 150 TeV). The 
ratio of /pp to the the acceleration length, using the pulsar wind density p determined by the 
value of /?, t a from Eqn (28) and magnetic field from Eqn (2), is 

^ =3.7 x (£) £) £) (29) 

Thus, the requirement /pp > /„ gives an upper limit on the proton energy E- r«v- Since it is 
expected that 7 » 1, possibly as large as 10® (Kennel and Coroniti 1984a), this limit will 
not be important in most cases. 

Protheroe (1983) has studied pair production and photopion production by > 10 12 eV 
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protons in UHE sources and gives attenuation times in blackbody radiation fields of 7 bb « 
10 4 — 10 s K from the stellar companion. The attenuation times for photopion production, 
important for protons with E > 10 15 eV, range from ss 10 j s(7bb = 8 x 10 4 K) to « 
10*8 (Tee = 10 4 K) and attenuation times for pair production vary from « 10 4 — 10® s for 
the same temperature range. With the shock velocity tii « c, as in the case of a relativistic 
pulsar wind and the magnetic field from Eqn (2), the strong shock acceleration time from 
Eq. (28), t a = 10 -6 s E^cv cm), is very short compared to the attenuation 
times above, even for energies > 10 18 eV. 


b) Maximum Acceleration Energy 


The maximum energy to which (charged) particles can be accelerated in the shock is de- 
termined by the balance of the energy gain rate with losses. The time t a needed for the 
shock to accelerate particles to an energy E was given in Eqn (28). The energy gains of 
the particles crossing the shock compete with energy losses through radiation and inelastic 
collisions and with diffusion away from the shock. For protons, diffusion is the dominant loss 
process since radiative losses are negligible even at the higher energies. Synchrotron losses 
are unimportant in limiting the maximum proton energy in nearly all cases (see below). The 
time to diffuse a distance r is t d = r 2 /D , where D is the diffusion coefficient and r is usually 
taken to be the shock radius. Equating the acceleration time, f„, to the diffusion timescale 
t d gives the maximum acceleration energy. With the minimum value of D, the maximum 
proton energy is 


j?m ax 


3((-l) 

N «{ + >) 


e£,r, 


(30) 


where the characteristic energy e B,r « 10 7 TeV in the case where r = r,. Requiring 

the particle gyroradius to be less than the shock radius results in the same limit, to within 
a factor of 2 or 3. The maximum energy E is approximately equal to the potential drop 
across the pulsar polar cap, <f> p = er 8 fi 2 Z? 0 /c. Because <j> p also equals the potential drop 
across the shock, the same E would result from drift acceleration in the perpendicular 
case. The electron acceleration will be limited by synchrotron losses, so their maximum 
energy is determined by equating the acceleration time in Eqn (28) to the synchrotron loss 
timescale. 
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1. Static atmosphere 


In the case of a shock in a static isothermal atmosphere, the proton energy will be limited by 
diffusion downstream into denser parts of the atmosphere where collisions become dominant. 
The scale length for diffusion is therefore the atmospheric scale height h rather than the 
shock radius r«r, » h. The condition determining the maximum acceleration energy is 
therefore t a = t j « h 7 /D. With the shock velocity ui = c, the magnetic field from Eqn (2) 
evaluated at r, and scale height from Eq. (6), the maximum proton energy is 


ET 


1.02 x 10® TeV 


1 


-r./«) s/ ‘ (31) 


The electron maximum energy is limited by synchrotron radiation losses whose timescale 
is t, = 200 sB~ 2 E^ly. Equating this timescale to the acceleration time t a and taking the 
magnetic field from Eqn (2), the maximum electron energy is 


= .033 TeV Bn' 2 P m (-) 1/2 (£-) ,/2 (1 - -) ,/2 

r. fiQ 


3(*-l) 


« V(e + i) 


(32) 


Figure 6 shows the regions of parameter space defined by Eqs. (13), (29) and (31 ) where a 
strong (£ = 4) pulsar wind shock in & stellar atmosphere could accelerate protons to energies 
above 10 and 10 4 TeV. The limits have been plotted for A//M© = r.//E© = nf K t — 1, 
although the dependence of the L & limits on these quantities is not very strong. There is 
a very restricted, but finite, parameter space (defined by the shaded strip) where protons 
can be accelerated above 10 TeV and no parameters for which acceleration to 10 4 TeV is 
possible. The general conclusion is that protons up to around 10 TeV can be produced by a 
pulsar wind shock which forms due to confinement by the atmosphere of the companion in 
a close binary. 


2. Stellar wind or magnetosphere 

In the case of confinement by a stellar wind or magnetosphere, the diffusion length scale 
which determines the maximum proton energy depends on how close the pulsar wind shock 
is to the companion star. If r, << a — r m , then the radius of the shock will be approximately 
r $ and = r 2 /D, but if r, ss a - r. with the shock near the companion star, then the shock 
radius will be approximately a - r, and tj = (a - r M ) 2 /D. The maximum proton energy in 


13 


these two cases is 


E™ = 1.2 x 10 7 TeVfi 12 P^ 2 , 


^ 4(4 + 1 ) 


f i, 


r, << a - r. 




(33) 

, % a - r. 

Note that in the case r, << a — r,, the r, dependence in ££“* (cf. Eqn [30]) drops 
out, because B t oc 1/r,, so that the maximum proton acceleration energy depends only 
on pulsar parameters. In fact, E™** cx LJ . The synchrotron loss timescale for protons, 
t , = 5xl0 ls sZ? _2 £f e v, is a factor (m p fm e ) 4 longer than for electrons. The maximum energy 
obtained by equating the synchrotron loss timescale to the acceleration time in Eqn (28) is 
lower than the maximum energy given in Eqn (33) when ( r s/P-_ )'^ 6 < 3.6. Thus, 

synchrotron losses are important in limiting the proton acceleration energy only in the case 
of millisecond pulsars with high magnetic fields. These pulsars, if they exist, would have 
spin-down timescales of 10 — 100 yr. 

As before, the maximum electron energy is determined by the balance of acceleration 
time and synchrotron loss time, but in the stellar wind case depends on the location of the 
shock: 


£T* = •033TeVBr2 ,/2 ^n S (r) ,/2 (| L ) ,/2 (-) ,/2 

/i© 


3 ( 4 - 1 ) 
\ 4(4 + 1 ) 


(34) 

r. /f© a 

The general conclusion here is that a pulsar wind shock due to confinement by a stellar 
wind (regardless of whether or not it is induced by the pulsar wind heating) or magnetosphere 
is capable of accelerating protons to energies of 10 TeV if P m ,B ,~ 2 ^ 2 < 900 and 10 4 TeV if 
Pts*.Bx 2 < 30. These limits translate, from Eq. (1), into pulsar wind power limits of 
Li> 6 x 10 31 erg/s and Lj > 5 x 10 37 erg/s. As mentioned in Section lib, induced wind 
solutions where M is large enough that the Alfven radius is inside the pulsar light cylinder 
would be unstable. These unstable solutions would occur for longer period pulsars with large 
light cylinder radii. If the accretion rate from the stellar mass loss is around the Eddington 
limit, then the requirement for a stable solution from Eqn (3) would be P na Bu^ < 31. This 
requirement would essentially restrict ihe parameter space to only those systems capable of 
accelerating protons to 10 4 TeV. 


IV. High Energy Gamma-Ray Emission 

Protons accelerated at the pulsar wind shock could produce 7-rays and neutrinos bv decay 
of pions, which would result when the protons interact with surrounding material (or con- 
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ceivably with ambient photons if their density is high enough). The protons that diffuse 
away from the shock are no longer confined in the pulsar wind and can cross the contact 
discontinuity to reach the companion star or other target material. Photons produced in 
interactions in the companion star wind or atmosphere would likely be modulated at the 
orbital period. If the shock is near the companion then the target subtends a relatively large 
solid angle at the proton source region around the pulsar wind shock (see Fig. 1). This 
model therefore may have a higher efficiency for converting accelerated proton luminosity 
into 7-rays than models in which protons are accelerated at the compact object. In the latter 
case the target subtends a much smaller solid angle at the proton source. 

The production spectrum of photons at the source is given by a convolution of the 
differential spectrum of accelerated protons, <f> p (E p ) = K E~ a , with the target density, the 
production cross section for 7r°’s and the ir° decay spectrum (Gaisser, 1988). Explicitly, the 
rate of photon production per target H atom is 

„ t E ?“ , „ 4E* 


dn, _ o [ E ?“ a v ( Bp a * 

dE y 2 k Ep Je, El 


<i> P {E p ). 


(35) 


Here x = E r /E p , and the momentum distribution of pions produced when protons collide 
in hydrogen is given by 


d?i p _ r o 


/(*: 


d E w E, 

The extra factor of 2/E r in Eqn (35) comes from the photon distribution in 7r° — ♦ 27. The 
accelerator is assumed to produce a power law spectrum of protons (differential index o) 
with a sharp cutoff at maximum energy E* 

One integration by parts leads to 
dn. 


r’m&x 
j p ■ 


dn 7 2 A j f l dz 1 


(36) 


where = E^/E\ 
of photons as 


max 
P * 


A similar manipulation gives an expression for the integral spectrum 

2 K 1 


tf,(> £,) = 




(37) 


where 


o(a-l) (£,)•"* " 

f b" + (o-lK)-M(*,r'l/(*)^- ( 36 ) 

JZ-y U 1 Z 

The function f{x) is determined by accelerator data. Berezinsky & Kudryavtsev (1989) 
give an explicit expression whi_u we use for numerical calculations: 

f(x) = 0.61 (1 - i) 3 - 5 + 0.46 exp(— 181). (39) 
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In the high energy limit (x^ -4 0), the function Z pr o(x y ), plotted in Figure 7, approaches 
the spectrum-weighted moment of the inclusive cross section for production of ir 0, s, 

Z pr o(0) = [ x°~ 2 f(z)dx. (40) 

The moments are given in Table 1 for a range of spectral index. 

The normalization constant for the spectrum of accelerated particles, which appears in 
Eqs. (35)-(37) is 


K = € p x Li 
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\Je^ E p 


SlY 1 
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(41) 


where t p is the efficiency for proton acceleration. This efficiency may be estimated by the 
solid angle, Q,, which the shock subtends at the pulsar (i.e. the solid angle over which the 
pulsar wind is confined): 




(42) 


which is only approximate because the shock is not spherical. 

The total yield of photons in the target region (before reabsorption) for a thin target is 


K,(> £,) = 


(43) 


where A is the proton interaction length and A A' the target thickness. For large pathlengths 
in the target region, A AT in Eq. (43) is replaced by the nucleon attenuation length, A = 
A/(l — Znn), which takes account of the fact that a fast nucleon interacts several times 
before losing all its energy. The quantity Znn depends on a, and is included in Table 1 . To 
estimate the magnitude of the photon flux seen by an observer at Earth, we next need to 
include a factor for the degree of beaming. This factor is l/(Afi</ 2 ), where Afi is the solid 
angle into which the accelerated protons are emitted and d is the distance to Earth. 

Finally, we must account for the duty cycle, A <j>^, which represents the fraction of the 
orbital period during which photons are emitted toward the observer. There are large un- 
certainties in estimating this factor. In particular, the signal that is emitted from the source 
depends crucially on the extent to which photons are reabsorbed in the target region. If the 
target is the companion star, we get the Vestrand & Eichler (1982, 1984) picture in which 
photons are only emitted along lines of sight through the limb of the companion. This sce- 
nario results in a rather low average signal, though broadened somewhat if the proton source 
is extended over the whole region of interface between the companion wind and the pulsar 
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wind (see Fig. 1). At the other extreme, one can imagine a “leaky source” in which charged 
particles, after acceleration, are confined by diffusion in turbulent magnetic fields until they 
have lost all their energy in collisions, but the line of sight thickness of this target region is 
thin enough so that all produced photons escape. In this case the beaming factor must be 
assumed to have its maximum value of AQ = 4 it. An intermediate case occurs if the target 
is the companion star wind, compressed by its interaction with the pulsar wind (see Fig. 1). 
This, however, would be a thin target case (A A' ~ 10 -3 gem -2 for PSR 1957+20), with a 
relatively small yield of photons. 

Assembling all the factors, the integral photon flux at Earth, averaged over the orbital 
period of the binary source, is 

AA* AX 

*•» (> £») = ^ (44) 


where 


2 z „ .(*,) ( dE p r 1 

(£p )-v ’ 


(45) 


^ a (a - 1) (£,)•-» [J Ea 

relates the photon spectrum to the parent spectrum that produces it. The factors in this 
equation have been grouped in this way for comparison with earlier estimates of the relation 
between the observed signal and luminosity at the source, especially the estimate of Hillas 
(1984) for Cygnus X-3. In addition, for > 10 H eV, absorption of photons in the microwave 
background due to 77 — > e + e~ must be accounted for. In the next section we use these ideas 
to estimate the signal from various sources. 


V. Binary Sources of VHE and UHE Emission 

There are known binary systems where acceleration in a pulsar wind shock could be occur- 
ring. The low mass X-ray binary Cygnus X-3 has been reported as a source of TeV and 
possibly PeV 7-rays at the 4.8 hr X-ray period (see Watson 1986 and Goodman 1989, for 
review). Both accretion and pulsar rotation have been suggested as the power source in this 
system, but the energy implied by the TeV and PeV 7-rays may favor (and perhaps require) 
a pulsar. There are also a number of binary systems known to contain spinning-down pul- 
sars. In particular, the recent discovery of an eclipsing radio pulsar (Fruchter et al. 1988) 
has generated interest in the interaction of a pulsar wind with a companion star. We wiL 
discuss the application of the pulsar wind shock acceleration model to these systems and 
give predicted fluxes of high-energy 7-rays. 
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Table 2 lists binary pulsars with observed values of period P and period derivative P, 
as well as some parameters for the Cyg X-3 system. Source distances (except for Cyg X-3) 
were determined from radio pulse dispersion measures assuming a mean interstellar electron 
density (n e ) = .03 cm -3 . The companion masses in the binary pulsar systems are better 
determined than the Cyg X-3 companion mass. The value 0.5Af© for the Cyg X-3 companion 
listed in Table 2 is the largest main sequence star which will fit inside its Roche Lobe 
(Patterson 1984), but the actual mass could be larger if it is a Helium star. Also listed in 
Table 2 is the pulsar luminosity incident on the companion star, L , = eLj. 

Since the predicted maximum proton acceleration energy in the pulsar wind shock model 
is proportional to the square root of the pulsar dipole luminosity, E™** may be determined 
for the case r, << 1 — r, directly from the spin-down luminosity (if it is known), without 
having to know the pulsar period and magnetic field strength separately. In cases where the 
shock is close to the companion, the maximum proton energy will be reduced from this value, 
depending on the value of the shock radius. Determination of the shock radius, however, 
requires a model for the pulsar wind confinement, several of which were discussed in Section 
II. Which confinement model or models to apply to the sources in Table 2 is not clear. None 
of the sources has L. > Lead required for an induced wind driven by radiation pressure (cf. 
Section Il.b.l), although this is not certain in the case of Cyg X-3 as the companion mass 
could be less than 0.5M©. We could therefore use the thermal wind model, which does not 
require L m > Lead- According to Eqn (21), the shock solutions in this model are determined 
only by the wind terminal velocity Voo/tW, which implicitly sets a value of M. Eqn (22) sets 
an upper limit on from the pressure balance requirement. There is also an upper limit 
on M from Eqn (3) for non- accreting sources like binary radio pulsars whose period and 
magnetic field strength are known. This upper limit on M is, strictly speaking, a limit on 
the neutron star accretion rate; the mass loss rate of the wind could be larger (but probably 
not much larger) if some of the mass escapes the system. The maximum values of A/ for 
the sources in Table 2 are listed in Table 3. We may therefore set a lower limit on from 
Eqns (3) and (19). This gives a range of possible shock radii from the thermal wind model 
for the sources in Table 2. Given the uncertainty in the shock radius in these sources, we list 
the maximum proton acceleration energy in Table 3 without the geometric reduction factor, 
and compute gamma-ray fluxes using these values, but we will consider the above limits in 
a few of the sources. 
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We make the following assumptions for the parameters defined in Section IV to tabulate 
the x°- decay photon flux from Eqn (44). Accelerated protons will escape downstream of the 
shock (*\e., away from the pulsar) and so have a preferred direction. However, they will not 
be highly beamed, since the shock is an extended structure. We therefore take the pn ton 
solid angle AQ = 1 sr. Within this solid angle, there are paths which encounter diff -.*ent 
thicknesses of target material, ranging from those going through the center of the companion 
star to those traversing less than 1 g cm -2 . The gamma-ray production efficiency increases 
with target thickness up to a maximum around 50 - 100 g cm -2 , and then decreases due 
to pair production attenuation. The peak gamma-ray signal will occur along those lines of 
sight where the gamma-ray production efficiency is maximum, where we can assume AX = A 
(thick target approximation). We also assume the proton acceleration efficiency, t p = | (its 
maximum value), and the gamma-ray duty cycle, A<t>-, = 1 (which effectively gives a peak 
flux). These are all n iltiplicative factors with which the predicted flux scales linearly. We 
show the results in Table 3 for two values of spectral index, a. Note that for steeper spectra 
the >100 MeV flux is larger, but the > 1 TeV flux is much lower than for the flat spectrum. 
This is a simple consequence of the fact that for steeper spectra more of the luminosity is 
dumped into the low energy portion of the spectrum. 

a) Cygnus X-3 

Some previous estimates of the power required to explain signals from Cyg X-3 have been 
in excess of 10 39 erg/s (Hillas 1984; Nagle, Gaisser and Protheroe 1988), and the required 
maximum proton energy is E™** > 10 16 eV. The fact that this luminosity exceeds the 
Eddington limit for steady accretion onto a neutron star (~ 10 M ergs -1 ), favors a model in 
which the power source is rotational energy release by a fast pulsar (Bignami et al. 1973, 
Vestrand and Eichler 1982). A period of 12.6 ms in TeV ->-rays from Cyg X-3 has been 
reported by the Durham group (Chadwick el al. 1985). The periodicity, not yet confirmed by 
other groups (Ramana-Murthy 1989), appears sporadically in the signal and measurements 
of changes in this period over 7 years give a period derivative Pm 3 x 10~ 14 s s -1 (Turver 
1989, private comm.). These values of P and P give a magnetic field strength of 6 x 10 n G 
and dipole luminosity 6 x lO^ergs -1 . 

It is interesting that the value of M < 9.3 x 10 ,8 gs -1 derived from Eqn (3) using the 
values of P and P reported by Chadwick et al. (1985) from TeV observations is just under 


19 


the upper limit, M < 10“ gs” 1 , on the Cyg X-3 system mass loss from the observed orbital 
period derivative (Bonnet-Bidaud and Van der Klis 1981). Since the shock location and 
maximum proton acceleration energy depend on the pulsar dipole luminosity, we do not 
need to know the pulsar period and magnetic field separately in applying the pulsar wind 
shock model to the Cyg X-3 system. We do need to know or assume values for the binary 
orbital period, companion mass and radius, separation and distance. The values we take for 
these parameters are shown in Table 2. 

The maximum proton acceleration energy, E^**, scales with L l J 7 . Assuming the dipole 
luminosity Ld = 6 x lO^rgs -1 implied by the values of P and P observed by Chadwick et 
al. (1985) gives E = 3 x 10 4 TeV. Using the formula in Eqn (44) for two values of proton 
spectral index a, the predicted flux of > TeV 7-rays from the source is 3.4 x 10~ 9 ph cm -2 s -1 
for q = 2 and 3.0 x 10 -1, ph cm -2 s -1 for o = 2.7. The predicted flux of > 100 MeV 7-rays is 
3.4 x 10 -5 ph cm -2 s -1 for a = 2 and 1.9 x 1 0“ 4 ph cm -2 s _1 for a = 2.7. In order to compare 
with phase-averaged observed fluxes, our peak flux predictions should be multiplied by the 
observed 4.8 hr duty cycle, ss 0.1. With the assumptions made, the a = 2 TeV flux 
predicted for Cygnus X-3 is greater than has been observed (Dowthwaite el al. 1984). The 
> 100 MeV fluxes tabulated are consistent with the SAS 2 measurement (SAS 2 is 1.1 x 10 -5 
above 35 MeV, Lamb et al. 1977, Fichtel et al. 1987 and earlier references therein)) ->ut 
are somewhat above the upper limit from the COS B observations (Hermsen et al., 1987), 
which is about 10 -6 for > 70 MeV photons. The primary reason that the predicted fluxes 
are high in the pulsar wind acceleration model is because of the much higher efficiency for 
generating a gamma-ray signal from Cyg X-3 than the Vestrand-Eichler model Since the 
accelerator is closer to the target material, a greater fraction of the accelerated protons can 
interact to produce gamma-rays. Consequently, the fraction of the pulsar spin-down power 
which is converted to accelerated protons need not be as large. Alternatively, emission from 
the source may be sporadic at the highest energies. 

Applying the thermal wind model to the Cyg X-3 system, we find that there is no range 
in wind terminal velocity which, as discussed above, gives a solution to Eqn (8) for the shock 
radius and satisfies the limit on M. Physically, the ram pressure of an induced thermal wind 
with M < A/„,*x (at least with the choices for parameter values made in Section II.B.2) 
is not large enough to stand off the pulsar wind above the companion surface. Either the 
parameter values chosen are not appropriate or the thermal wind model of Cheng (1989) does 
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not work for this system, possibly because it is close to the Eddington limit. Alternatively, 
it is possible that only part of the pulsar wind exerts pressure on the companion star wind 
(Rasio, Shapiro and Teukolsky 1989). In any case, because Cyg X-3 is a close binary the 
shock radius »vould be of the order of r,/a « 0.5 and the geometrical reduction factor in 
Eqn (33) for ££“* is close to unity. 

b) Binary Pulsars 

We have calculated the maximum proton energy in the pulsar wind acceleration model for 
nine known binary radio pulsars, listed in Table 2 with their period P, period derivative 
P, orbital period companion mass M c and distance d (cf. Dewey et al. 1986). Table 3 
gives the values of A/ max and E™**, as well as the predicted ',-ray fluxes from Eqn (44) for 
o = 2 and o = 2.7. Several sources have predicted fluxes of > 100 MeV 7-rays which are 
above the sensitivity threshold of 5 x 10 -8 phcm _2 s -1 of the EGRET detector (Kanbach et 
al. 1988) on the Gamma-Ray Observatory (GRO). Because of the dependence cf maximum 
proton energy on pulsar luminosity, these are also the systems capable of producing 7-rays 
above 1 TeV. One of these, PSR1957+20, is the eclipsing millisecond pulsar which may be 
evaporating its companion by means of an induced wind (Phinney et al. 1988, Kluzniak 
et al. 1988, Cheng 1989). A mass loss rate of M ss 10 16 gs -1 is suflcient to evaporate it 
completely in the pulsar’s lifetime. Another system, PSR1855+09, contains a 5 ms pulsar 
and is relatively nearby, giving a predicted flux > 100 MeV well above EGRET sensitivity 
(for the assumptions we have made). 

Since the r, solutions in the thermal wind model are parameterized by the wind terminal 
velocity, we must specify a value of v,*, in the solution for r,. Since r./a << 1 in all of these 
systems, the upper limit on v ^ from Eqn (22) is around 0.29 In several systems, notably 
PSR1913+16, the companion is probably another neutron star, which would not be expected 
to have a wind. In that case, the magnetic pressure of the neutron star companion or ram 
pressure of another pulsar wind would balance the pulsar wind, with the shock location 
determined by Eqn (23). 

In the case of the eclipsing pulsar PSlvi 957+20, it is not clear whether a pulsar induced 
thermal wind is capable of supporting the pulsar wind above the stellar surface. The mass 
loss rate in the X-ray heated thermal wind model, using Eqn (17), is Af ~ 6 x 10 I3 gs -1 
(Emmering and London 1989), assuming Uqo = t>e*c- The solution for the shock location 


21 


from Eqn (18) gives r t /a ~ .985, which is inside the star. In the model of Cheng (1989), a 
value (voo/tw) = 0.1 puts the shock at the observed eclipse radius, r,/a ss 0.6. The formula 
in Eqn (10) gives 1 — r,/a ss 9.3(r„/a) for (voo/*>«c) = 0.1 when A w >> 1, so in this case the 
shock stands off about 9.3 stellar radii firnn the companion. This value of r,ja reduces 
only slightly from the value in Table 3 to 250 TeV, which will not change the > 1 TeV flux 
estimate. However, (v<»/tw) = 0.1 gives a value A/ = 6x 10 17 gs _l from Eqn (19), which 
is much larger than the rate (given ir. . ble 3) at which the pulsar would begin to accrete 
mass. Increasing (««, /»<*c) to 0.5 woe' ' catisfy the limit on Af, but the induced wind is then 
barely able to stand off the pulsar v, i.,d above the companion surface. In this case, JE£“ 
would be reduced to 20 TeV, with a resulting decrease in the > 1 TeV 7-ray flux by a factor 
of 3 for q = 2. 

In the case of PSR185O+09, a value of (t7ooA«c) > .057 is required to keep the induced 
M in the thermal wind model below A/ m , v . However, we require (v^Aw) < .29 for pressure 
balance with the pulsar wind, so there is a range of allowed solutions. Even in the most 
favorable case, M = Af miut , the solution for the shock radius is then r,/a = 0.95 and is 
reduced to only 4 TeV, lowering the predicted flux to 0.1% of the value in Table 3. 

VI. Discussion 

We have proposed a model for particle acceleration in VHE and UHE 7-ray sources which 
contain rapidly spinning, non-accreting pulsars. Considering the physics of the first-order 
Fermi mechanism at the pulsar wind shock, acceleration to energies above 10 17 eV is in 
principle possible (e.g. if the source contains a 10 ms pulsar with a 10 12 G magnetic field). 
Furthermore, the power source is the pulsar spin-down, which could supply superEddington 
luminosities in accelerated particles. Models for UHE sources in which the power source is 
accretion have a luminosity limit of L < 2 x lO^ergs'^Ro/R), where R is the radius at 
which the accretion energy is converted particle energy. These models also have difficulty 
accelerating particles to the energies needed to produce UHE 7-rays (see Harding 1989 for 
review). The pulsar wind model is able to generate the combination of UHE 7-ray emission 
and possibly superEddington luminosity of Cyg X-3. In addition, the accelerator-target 
configuration in this model may be more favorable, so that the power required to produce 
a given signal is not so severe. This acceleration model would not apply to sources like Her 
X-l and Vela X-l, which have been detected at VHE and UHE energies, but are known to 
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be accreting neutron stars. In these cases though, the power available through accretion is 
sufficient to account for the observed 7-ray fluxes. 

If confined pulsar winds result in acceleration of particles, then some binary pulsars which 
nave been detected by their radio pulses could be 7-ray sources. Recently, emission at TeV 
energies has been reported by von Ballmoos et al. (19S9) and de Jager et al. (1989) from 
two binary pulsars. Data from PSRl 957+20, folded with the 9 hr orbital period, shows a 
peak in the phase plot at the position of the LA Lagrange point. The reported flux above 2 
TeV is 1.1 x 10~ 9 photons cm -2 s -1 , much higher than our predicted fluxes in Table 3. This 
difference is due to the small 7-ray solid angle of Aft ss .009 derived by von Ballmoos et 
al. (1989) from the width of the phase peak (A^ = .018), whereas we have taken Aft = 1 
sr to compute the fluxes in Table 3. A peak at the LA position in the orbital phase plot 
of PSRl 957+20 may also have been seen at > 100 MeV in COS-B data (von Ballmoos et 
al. 1989). A TeV signal was also reported by de Jager et al. (1989) at the 5.4 ms pulsar 
period of PSR1855+09 at a marginal significance level. Signals which are observed to be 
pulsed at the pulsar period, however, must originate from acceleration near or within the 
pulsar magnetosphere and would not be expected from particles accelerated at the pulsar 
wind shock. 

The currently favored mode) for the origin of binary systems containing short period 
pulsars is spin up by an accretion disk (Alpar et al. 1982). The accretion spin-up period 
depends on the neutron star magnetic field in this model, with a spin-up period in the 10 ms 
range requiring a field around 10 9 Gauss. This evolution model accounts quite well for the 
observed binary pulsars, which all have low magnetic fields. The high luminosity implied 
by the observed Cyg X-3 gamma-ray flux requires not only a short period but also a high 
magnetic field, incompatible with an accretion spin-up evolution. A fast pulsar as the power 
source in Cyg X-3 would therefore have originated in a relatively recent supernova explosion 
that occurred within the pulsar spin-down time of 7000 yr. Confined pulsar wind sources of 
this type would be more rare than the accretion spun-up binary pulsars, which have ages of 
around 10® yr. 

In this paper, we have examined pressure balance only along the line between the pulsar 
and the companion. In other directions, the balance point, and thus the shock location, 
would be further from the companion. The details of the particle acceleration and gamma- 
ray production (predicted light curves, etc.) will require a full solution to the pressure balance 
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and the sh&p. of the pulsar wind shock. In addition, such a solution may allow asymmetric 
pulsar-induced winds in cases where there is no pressure balance above the stellar surface 
Thus, a thermal wind ram pressure too low to support the pulsar wind above the surface 
facing the pulsar does not altogether prevent ablation from the limb of the star. 

We also have not addressed the question of the source of protons which are accelerated at 
the shock. One possible source of protons is the pulsar wind, which in most models contains 
only electrons and positrons, but could conceivably carry some ions. Another possibility 
is that some of the material from the companion star wind or atmosphere may mix into 
the pulsar wind just downstream of the shock, through Rayieigh-Tayior instabilities at the 
contact discontinuity. Such a mechanism for mixing has been investigated in the context of a 
pulsar wind confined by a supernova remnant (Gaisser et at. 1989, Harding et at. 1989). In 
the binary system, though, the distance of the pulsar wind shock to the contact discontinuity 
is not as easily determined, since it would require a steady-state solution for the flow of the 
pulsar wind around the companion. 
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Table 1 

Parameters for Photon Production 


a 

2.0 

2.2 

HI 

2.6 

2.8 

Znn 

m 


0.34 

0.31 

0.29 


0.17 


0.066 

0.048 

0.036 

(includes 6% 17) 







Table 2 

Binary Pulsar Parameters 


PSR P Log(P) P b M c a d L d L m 

(ms) (days) (A/©) (cm) (kpc) (erg/s) (erg/s) 


1913+16 

59 

-17.1 

0655+64 

195.6 

-18.2 

1831-00 

520.9 

-17 

1855+09 

5.4 

-19.8 

2303+46 

1066.4 

-15.4 

1953+29 

6.1 

-19.5 

0820+02 

864.9 

-16 

1957+20 

1.61 

-19.9 

1620- 26 

11.1 

-18.1 

CYG X-3 

12.6 

-13.5 


0.32 

1.4 

1.9(11) 

1.03 

1 

4.0 (11) 

1.81 

0.1 

5.0 (11) 

12.33 

0.3 

1.9 (12) 

12.34 

2 

2.4 (12) 

117.35 

0.3 

8.4 (12) 

1232.4 

0.3 

4.0 (13) 

0.381 

0.02 

1.7(11) 

191.44 

0.35 

1.2 (13) 

0.2 

0.5 

1.2(11) 


4.33 

1.58 (33) 

7.1 (28) 

0.27 

3.45 (30) 

4.5 (25) 

3.13 

2.90 (30) 

1.1 (26) 

0.44 

4.12 (33) 

5.6 (27) 

2.0 

1.34 (31) 

3.2 (24) 

2.92 

5.71 (33) 

3.8 (26) 

0.79 

6.33 (30) 

1.8 (22) 

1.0 

1.24 (35) 

1.2 (32) 

2.1 

2.38 (34) 

7.3 (26) 

10 

6 (38) 

2.1 (37) 
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Table 3 

Gamma- Ray Flux from Binary Pulsars 



PSR 

Mmmx 

(g /•) 

r./a 

(TeV) 

#,(> ITeV) 
(ph cm -2 s -1 ) 
a = 2.0 

t'; 

CM 

II 

0 

*.,(> lOOMeV) 

(ph cm -2 s' 1 ) 
o = 2.0 a = 2.7 

1913+16 

2.0 (13) 

1.34 (-2) 

50.3 

4.5 (-14) 

2.5 (-16) 

7.3 (-10) 

2.6 (-9) 

0655+64 

2.1 (11) 

7.25 (-3) 

2.35 

— 

— 

5.9 (-10) 

1.6 (-9) 

1831-00 

2.8 (11) 

1.25 (-2) 

2.15 

— 

— 

3.6 (-12) 

9.2 (-12) 

1855+09 

4.1 (13) 

2.32 (-3) 

81.2 

1.2 (-11) 

1.0 (-13) 

1.6 (-7) 

6.8 (-7) 

2303+46 

1.9 (12) 

9.78 (-4) 

4.64 

— 

— 

3.7 (-11) 

1.0 (-10) 

1953+29 

6.0 (13) 

5.17 (-4) 

95.6 

4.2 (-13) 

3.3 (-15) 

5.5 (-9) 

2.2 (-8) 

0820+02 

7.9 (11) 

1.08 (-4) 

3.18 

— 

— 

1.2 (-10) 

3.2 (-10) 

1957+20 

6.8 (14) 

6.17 (-2) 

445 

8.9 (-11) 

6.2 (-13) 

8.9 (-7) 

3.9 (-6) 

1620-26 

3.4 (14) 

3.51 (-4) 

195 

3.6 (-12) 

2.7 (-14) 

4.2 (-8) 

1.7 (-7) 
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FIGURE CAPTIONS 


Figure 1 - Schematic view of pulsar wind * ock formation in a binary system for the case of 
confinement by the companion star wind. 

Figure 2 • Distance of the pulsar wind shock from the pulsar vs. companion radius for 
the general case of confinement by a companion star wind. Curves are labeled 
with values of the dimensionless parameter, A w , defined in Eqn (11). The curved 
right-hand boundary of the solutions is where r, = a - r m . 

Figure 3 - Distance of the pulsar wind shock from the pulsar vs. companion radius for 
the case of confinement by a companion star wind driven by radiation pressure. 
Curves are labeled with values of 

Figure 4 - Distance of the pulsar wind shock from the pulsar vs. companion radius for the 
case of confinement by a companion star thermal wind. Curves are labeled with 
values of v^/t 

Figure 5 - Distance of the pulsar wind shock from the pulsar vs. companion radius for the 
case of confinement by a companion star magnetosphere. Curves are labeled with 
values of the dimensionless parameter, log 10 (A m ), defined in Eqn (24). 

Figure 6 - Limits on = Lj/lO^ergs -1 and companion radius to binary separation ratio 
required for proton acceleration to energies above 10 and 10 4 TeV in a pulsar 
wind shock, for the case of pressure balance in the stellar atmosphere. 

Figure 7 - Relative efficiency of 7r°-decay photon production for power law proton spectra 
with index a = 2.0, 2.2, 2.4,?. 6, 2.8. The function Z pK o(x^) is defined in Eqns (38) 
and (40). The quantity = E y / E™**. 
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